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摘 要 :考虑 轴承 故障 动力 学 中 的 高 频 共 振 以 及 故障 激励 特点 ,提出 一 种 模拟 系统 高 频 共 振 和 过 渡 
式 故 障 激励 的 轴承 故障 动力 学 建 模 方法 。 基 于 赫兹 接触 理论 ,将 轴承 部 件 之 间 的 连接 简化 为 弹簧 - 
阻尼 结构 ,并 加 入 谐振 器 以 模拟 故障 的 宽频 激励 所 引起 的 系统 高 频 共振 ,从 而 建立 轴承 六 自由 度 动 
力学 模型 。 另 外 ,传统 的 故障 瞬 变 激励 方法 难以 反映 滚动 体 经 过 滚 道 故 障 的 过 程 ,考虑 到 滚动 体 与 
故障 接触 的 实际 情况 ,给 出 一 种 过 渡 式 的 故障 激励 方法 。 以 传动 链 平台 上 的 6304 深 沟 球 轴承 为 
生 * 例 ,对 轴承 外 圈 故 障 下 的 振动 响应 特性 进行 仿真 和 数据 处 理 ,其 仿真 结果 与 轴承 故障 实测 信号 对 比 
-相符 ,仿真 信号 故障 一 阶 特征 阶 次 与 理论 值 误差 约 为 0.0167% ,表明 模型 具有 一 定 的 准确 性 和 可 
OD 千 性 。 
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Abstract: Considering the characteristics of high frequency resonance and fault excitation in bearing fault 
dynamics ,a bearing fault dynamics modeling method is proposed to simulate high frequency resonance and 
transient fault excitation of the system. Based on Hertz contact theory ,the connection between bearing com- 
ponents is simplified to spring-damper structure ,and resonance changer is added to simulate the high-fre- 
quency resonance of the system caused by wide-band excitation of faults,thus establishing a 6-DOF dy- 
namic model of bearings. In addition ,the traditional fault transient excitation method can not reflect the 
process of rolling element passing through the raceway fault. Considering the actual contact condition be- 
tween rolling element and fault,a transitional fault excitation method is given. Taking 6304 deep groove 
ball bearing on drive chain platform as an example ,the vibration response characteristics under bearing 
outer ring fault are simulated and data processed. The simulated results are consistent with the measured 


signal pairs of bearing faults. The error between the first-order characteristic order of the simulated signal 
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and the theoretical value is about 0.0167% ,which shows that the model has certain accuracy and reliabil- 


ity. 
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齿轮 箱 是 传动 链 中 的 重要 结构 ,而 轴承 是 确保 
齿轮 箱 正常 运行 的 关键 组 成 部 分 。 对 于 深 动 轴承 来 
说 ,建立 有 效 的 轴承 故障 仿真 工程 模型 ,对 于 分 析 轴 
承 故障 成 因 有 着 重要 意义 ,为 轴承 故障 诊断 的 深入 
研究 提供 了 理论 基础 。 

在 轴承 动力 学 理论 研究 方面 ,诸如 疲劳 .刚度 、 
赫 效 接触 . 游 隙 影响 和 冲击 力 影响 等 ,选择 合适 的 研 
究 方法 是 轴承 高 效 建 模 的 重要 手段 。Sassi 等 通 
过 分 析 轴 承 的 动力 学 行为 ,将 轴承 建 模 为 内 圈 - 滚 动 
体 -处 圈 的 三 自由 度 模型 ,将 油膜 等 效 为 弹簧 阻尼 结 
构 动 ! 和 随机 噪声 并 考虑 经 过 故障 的 冲击 力 实 现 对 
族 限 铀 承 的 仿真 。Feng 等 "基于 赫 效 接触 建立 了 
畏 这 四 自由 度 模型 ,研究 了 轴承 内 外 圈 故 障 下 的 不 
周 困 应。Sawalhi 等 中 将 从 轮 箱 轴 承建 模 为 五 自由 
魔 模型 ,并 对 轴承 矩形 缺陷 进行 了 优化 。Patel 等 
将 滚动 体质 量 纳入 模型 建立 了 六 自由 度 轴 -滚动 体 - 
销 电 座 模型 ,轴承 部 件 之 间 为 赫兹 接触 ,用 矩形 缺陷 
的 局 寸 参数 计算 故障 激励 。Rehab 等 中 建立 了 轴承 
噬 骨 由 度 模型 ,研究 了 故障 尺寸 的 变化 对 故障 特征 
频率 处 幅 值 的 影响 。Mishra 等 "对比 了 轴承 五 自由 
麻 异 型 ,多 体 动力 学 模型 和 ADAMS 仿真 模型 的 故 
障 山 应 ,并 通过 实验 数据 证 明了 模型 的 准确 性 。 
sh 字 等 "利用 轴承 六 自由 度 模型 研究 了 干燥 与 渔 
浓 僵 况 下 轴承 故障 特征 的 区 别 。 李 志 农 等 "改进 
故障 模型 为 渐进 模型 , 利用 二 自由 度 模 型 仿真 验证 
了 改进 模型 的 正确 性 。 王 晨 等 29 研究 了 非 规则 故 
障 模型 及 轴承 故障 参数 对 系统 振动 响应 的 影响 。 罗 
茂 林 等 加 将 冲击 力 引 入 轴承 动力 学 模型 ,并 对 轴承 
故障 信号 的 双 冲 击 特征 进行 了 仿真 。 田 晶 等 52 利 
用 轴承 四 自由 度 模型 对 中 介 轴 承 多 点 故障 进行 仿 
真 ,并 对 比 了 缺陷 尺寸 .载荷 和 转速 比 对 故障 特征 的 
影响 。 

在 齿轮 箱 系 统 中 ,如 果 齿 轮 发 生 缺 齿 、 断 齿 等 故 
障 , 故 障 信号 主要 集中 在 低频 频段 ;而 轴承 故障 往往 
会 从 轴承 内 部 激发 周期 性 冲击 ,这 种 周期 性 的 冲击 
能 够 在 整个 系统 的 共振 频带 附近 产生 高 频 共振 ,所 
以 根据 这 一 特点 ,在 齿轮 箱 的 故障 诊断 中 利用 带 通 
滤波 器 选择 合理 频段 进行 分 析 , 能 够 获得 更 加 明显 


的 轴承 故障 信号 ,而 传统 的 轴承 二 自由 度 、 四 自由 度 
模型 难以 模拟 轴承 故障 这 一 特点 。 因 此 ,本 研究 在 
轴承 四 自由 度 动力 学 模型 ”基础 上 ,水 平和 竖 直 方 
向 上 额外 添加 了 两 个 谐振 器 来 激励 轴承 在 高 频 范 围 
内 的 共振 频带 ,并 且 给 出 一 种 更 加 合理 的 过 渡 式 故 
障 激励 方法 ,利用 数值 方法 求解 得 到 外 圈 故 障 的 轴 
承 振 动 啊 应 ,与 实测 轴承 故障 信号 进行 对 比 ,验证 了 
模型 的 准确 性 ,为 滚动 轴承 的 故障 诊断 提供 参考 。 


1 轴承 故障 动力 学 建 模 


1.1 考虑 谐振 器 的 轴承 集中 参数 模型 


常见 的 滚动 轴承 几何 参数 如 图 1 所 示 。 图 中 ， 
轴承 外 圈 直 径 为 D, ,内 图 直径 为 Di ;轴承 外 滚 道 直 
径 为 由 ,内 滚 道 直径 为 方 ;轴承 节 圆 直径 为 忆 , ,滚动 
体 直 径 为 D, ,轴承 径 向 游 际 为 C,。 


图 1 轴承 几何 参数 

Fig. 1 Geometric parameters of bearing 
集中 参数 模型 主要 考虑 轴承 部 件 的 平面 运动 ， 
轴承 滚动 体 和 深 道 之 间 的 接触 以 及 轴承 座 的 作用 简 
化 为 弹簧 阻尼 结构 ,轴承 部 件 围绕 轴 的 旋转 被 假设 
为 纯 滚动 ,忽略 了 摩擦 的 作用 ,六 自由 度 模型 如 图 2 
所 示 。 

模型 中 轴承 受到 的 轴 疝 力 可 以 忽略 不 计 , 重 点 
考虑 滚动 轴承 横 截 面 所 在 平面 内 的 两 个 相互 垂直 的 
径 向 振动 ,轴承 内 圈 固 定 在 刚性 转轴 上 ,轴承 外 圈 固 
定 不 发 生 转动 。 轴 承 的 自由 度 完全 集中 在 滚 道上 ， 
便于 确定 作用 在 内 圈 上 的 时 变 轴承 力 ,作用 在 轴承 
内 外 圈 上 的 力 取决 于 各 个 滚动 体 受 力 之 和 。 
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图 2 轴承 六 自由 度 模型 


T Fig.2 Six degree of freedom model of bearing 

司 设 定 径 向 平面 竖 直 向 下 为 了 轴 正 方向 ,水 平 向 
存 鸡 工 轴 正方 向 。 考 虑 到 滚动 体 与 滚 道 的 非 线性 接 
角 > 轴 承 六 自由 度 动力 学 模型 的 动力 学 方程 组 为 


@ mxX; tex thx; + Far =F,sin( wi) (1) 


EA my;+cey; +t+hky+Fyy=P+Fcos(wt) (2) 


moX, + (Co ten) Xo t+ (hot hy )%, — 


kyn 一 caya = 人 ar (3) 
my ,s+(cs, teo)y, + (kk, +hs)y, — 
poye 一 coya = Fhy (4) 


ma Xu toca(Xs -KX) +h (x -x,)=0 (5) 
LO =0 (6) 
式 中 :m, 和 m, 分 别 是 轴承 内 图 与 轴 的 等 效 质量 和 
外 团 与 轴承 座 的 等 效 质量 ;hk; 和 c; 是 轴 的 等 效 连接 
刚度 和 阻尼 ;6,, 和 分 别 是 轴承 座 在 水 平方 向 上 的 
等 效 连 接 刚 度 和 阻尼 ;k, 和 c 分 别 是 轴承 座 在 竖 直 
方向 上 的 等 效 连接 刚度 和 阻尼 ;x 和 y; 分 别 是 内 圈 
在 水 平和 坚 直 方向 上 的 位 移 ;x, 和 y, 分 别 是 外 圈 在 
水 平和 坚 直方 向 上 的 位 移 ;w 是 内 圈 转 动 的 角速度 
(r/min) Pe 是 作用 在 轴 上 的 负载 ;mv 和 Fy 分 别 是 
内 外 圈 与 滚动 体 接触 而 产生 的 弹性 恢复 力 在 水 平和 
坚 直方 向 上 的 分 力 ; ,是 旋转 轴 偏 心 引起 的 离心 
力 ;muyka ycn 分 别 是 水 平方 向 上 谐振 器 的 质量 、 刚 
度 和 阻尼 ;ms ,ks ,co 分别 是 竖 直 方向 上 谐振 器 的 质 
量 、 刚 度 和 阻尼 ;x,, 是 水 平方 向 上 谐振 器 的 水 平 位 移 
响应 ;ys 是 竖 直 方向 上 谐振 器 的 竖 直 位 移 响应 。 

设 第 一 个 滚动 体 的 初始 角 位 置 9, 为 0, 滚动 体 
个 数 为 N, 保 持 架 转动 频率 为 w。。, 则 第 ;个 滚动 体 
处 在 上 时刻 的 角 位 置 9 为 
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0 = to (7) 


内 滚 道 半径 为 R;, 外 深 道 半径 为 RR,, 则 保持 架 

wo = (8) 

轴承 在 运行 时 ,第 i 个 深 动 体 与 内 外 圈 接 触 产生 

的 接触 变形 y 取决 于 轴承 内 外 圈 的 相对 位 移 以 及 径 

向 游 际 C, ,内 外 圈 的 相对 位 移 可 由 几何 关系 确定 , 故 
轴承 滚动 体 接触 变形 量 y 的 取 值 由 式 (9) 决 定 , 即 

yi= (Xi—%,)sing, +(y;—Y,)cos0,-C, (9) 

考虑 到 实际 深 动 体 与 深 道 接触 变形 量 不 存在 为 

负 值 的 情况 , 即 深 动 体 与 内 外 深 道 接触 产生 的 接触 

变形 y; 恒 大 于 等 于 0, 如 式 (10) 所 示 。 


Yi, Yi >0 
» =| (10) 
0 ? Yi <0 
轴 安 装 误差 引起 偏心 产生 的 离心 力 为 
下 =em,w’ (11) 


其 中 e 是 偏心 距 。 

当 深 动 体 受 载 与 深 道 产生 接触 形变 时 ,单个 深 
动 体会 给 总 的 赫兹 接触 力 一 个 增 基 ,基于 赫兹 接触 
理论 得 到 轴承 内 外 深 道 的 等 效 接触 刚度 , 深 动 体 
的 赫兹 接触 力 在 X 轴 和 Y 轴 的 分 量 分 别 为 


N 
Pa = KY y,1. 5sing, (12) 
i=1 


N 
Fiy = 天》， y; 1. Scos0， (13) 
i=1 


根据 赫兹 接触 力 的 定义 可 以 发 现 , 深 动 体 被 等 
效 为 非 线 性 弹簧 ,通过 计算 时 变 的 赫兹 接触 力 可 以 
分 析 正 常 状态 下 轴承 运行 时 的 变 柔 度 振动 。 


1.2 ”过渡 缺陷 的 模拟 方法 


传统 的 故障 模型 为 瞬 变 激励 的 矩形 缺陷 , 当 滚 
动 体 进入 故障 的 瞬间 就 释放 出 所 有 形变 量 ,如 图 3 
(a) 所 示 ; 而 实际 的 轴承 故障 缺陷 形状 呈 半 椭 球 形 
凹 坑 , 当 滚动 体 进入 缺陷 时 ,形变 量 释 放 的 过 程 应 该 
是 缓慢 过 渡 的 '” 。 本 研究 的 故障 模拟 采用 一 种 新 
的 过 渡 故 障 激励 ,如 图 3(b) 所 示 。 

以 工程 中 轴承 常见 的 外 圈 故 障 为 例 , 基于 滚动 
体 与 缺陷 实际 接触 的 几何 关系 对 轴承 故障 进行 建 
模 。 将 故障 模拟 为 具有 椭 球 形状 的 四 陷 , 本质 上 是 
通过 游 陈 的 变化 来 表征 故障 , 建 模 时 能 够 在 仿真 参 
数 上 更 直观 地 反映 故障 的 不 同 程度 。 下 面 分 析 滚动 
轴承 滚 道上 的 局 部 故障 。 


(b) 过 渡 缺 陷 


图 3 故障 模型 示意 图 


Fig3 Schematic diagram of fault model 


滚动 轴承 在 运行 时 ,滚动 体能 和 否 接触 到 故障 底 


部 ,对 滚动 体形 变 的 释放 量 有 着 很 大 影响 ,滚动 体 进 


入 外 圈 故 障 处 局 部 放大 如 图 4 所 示 。 


R 
外 
图 4 滚动 体 与 外 圈 故 障 接触 示意 图 


“SFig.4 Schematic diagram of faulty contact between 


Ee 
《9? rolling element and outer race 


H, | cos’ (0,— pp,) -cos’( 


2 )] 


ChinaXiv 合 作 期 刊 485 


当 滚 动 体 进入 故障 区 域 时 , 深 动 体 与 深 道 接触 
瞬时 消失 , 当 其 离开 故障 区 域 时 ,滚动 体 与 滚 道 接 触 
恢复 ,接触 的 消失 和 恢复 来 源 于 局 部 故障 造成 的 变 
化 。 具 体 地 , 当 引 入 故障 后 ,轴承 滚动 体 接触 变形 量 
yi 的 取 值 将 由 式 (14) 决 定 。 

y;= (xX;—%,)sing, +(y;—Y,)cos0,-C.-H (14) 
其 中 右 是 轴承 滚动 体 在 经 过 故障 时 陷入 故障 的 距 
离 。 由 于 轴承 在 运行 时 的 内 外 圈 的 运动 不 同 ,对 于 
不 同类 型 的 轴承 故障 五 也 有 着 不 同 取 值 。 

深 动 体 进 入 故障 时 释放 的 形变 量 有 的 取 值 与 
深 动 体 与 故障 的 几何 关系 相关 。 考 察 故障 处 的 几何 
关系 , 当 故 障 深度 hh 足够 大 时 ,HH 的 值 为 滚动 体形 
变量 CG, 与 外 圈 形 变量 C。, 相 减 ; 当 故 障 深 度 h 较 小 
时 , 的 值 为 故障 深度 有 hh, 如 式 (17) 所 示 。 


rh 一 一 (15) 


-二 (16) 
8r,h=4h + 
0 二 可 (17) 
{i 8rih <4h + 
Via 为 轴承 故障 所 对 应 的 轴承 圆心 角度 ,可 根据 
几何 关系 由 式 计 算 。 


2L 
Wu TD (18) 


故 考虑 滚动 体 进 入 和 驶 出 缺陷 时 形变 量 的 释放 
过 程 , 当 轴 承 外 圈 发 生 故 障 时 ,滚动 体 陷 和 人 故 隐 的 距 
离 且 根据 式 (19) 取 值 。 


Bs 


H = 2/ Vou 
sin (>) 
0， 

式 中 :pu 为 外 圈 故 障 位 置 角度 ,由 于 轴承 外 圈 固 定 ， 
所 以 gi 为 一 常数 。 

将 故障 引入 轴承 六 自由 度 动力 学 方程 组 实现 轴 
承 故 障 动力 学 建 模 ,利用 四 阶 龙 格 库 塔 法 求解 可 以 
得 到 轴承 在 故障 状态 下 运行 的 振动 响应 。 


2 仿真 计算 与 实验 验证 


为 验证 轴承 故障 仿真 模型 的 正确 性 ,采用 实验 


8 他 = 三 c out 
cto -pn) on[ 加 i 


others 
室 传动 链 平台 的 齿轮 箱 实测 故障 信号 与 仿真 信号 进 
行 对 比 。 传 动 链 平台 由 驱动 电机 、 锥 齿轮 箱 系统 、 平 
行 轴 齿 轮 箱 系统 和 行星 齿轮 箱 系 统 组 成 ,故障 试验 
件 为 NSK 公司 生产 的 6304 轴承 ,如 图 5 所 示 。 

采用 光电 传感器 测量 轴 速 , 压 电 振动 加 速度 传 
感 器 测量 振动 信号 ,采样 频率 51. 2 kHz。 试验 件 的 
基本 参数 如 表 1 所 示 。 

实验 台 参 数 如 下 :内 圈 与 齿轮 轴 轴 质量 mm 为 


0. 9 kg, 外 圈 与 轴承 座 质 量 mm。 为 5.5 kg, 故 障 的 宽度 


486 


L 与 深度 用 均 为 1 mm。 实 验 中 齿轮 箱 没有 连接 负 
载 , 故 负 载 P 大 小 等 于 内 圈 与 从 轮轴 的 重力 。6304 
轴承 所 在 轴 的 转速 为 1 179. 31 r/min, 即 转 频 了/ 
为 19. 66 Hz。 


园 惟 占 轮 箱 | 


: 平行 


SH 齿轮 箱 齿轮 箱 
* 驱动 
电机 


PN 


图 5 传动 链 实验 平台 
Fig.5 Transmission chain experiment platform 
表 1 6304 轴承 试验 件 基本 参数 
Tab.1 Basic parameters of 6304 bearing 


0409v1 


4 参数 数值 
2 Di /mm 20 
Dy /mm 52 
Di/mm 9.525 

~ D,/mm 36 
SC 深 动 休 数 上 7 


5、 本 a 但 \ 了、 一 He 
忆 为 提高 轴承 模型 的 准确 性 ,通过 有 限 元 建 模 得 


到 实验 台 的 部 分 参数 。 利 用 ANSYS Workbench 
19e9 的 静 力 学 模块 对 轴 和 齿轮 箱 进行 有 限 元 分 析 
以 得 到 建 模 参数 ,如 图 6、 图 7 所 示 。 有 限 元 计算 得 
到 轴 的 刚度 态 为 3.419 x10 N/m, 齿 轮 箱 水 平方 向 
上 的 刚度 到, 为 3.677 x10”N/m, 竖 直方 向 上 的 刚度 
,为 9.326 x10 N/m。 


A.6304shaft 

Total Deformation 

Type:Total Deformation 

Unit:m 

Time:1 

2022/2/22 15:27 
0.00014624 max _ 
0.00012999 蝴 
0.00011374 
9.7492e-5 
8,1243e-5 
6.4995e-5 
4.8746e-5 
3.2497e-5 
1.6249e-5 Oe (my 
0 min 0.025 0.075 


到 6 轴 位 移 云 


Fig.6 Deformation nephogram of shaft 


器 
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C:box_shuiping 

Total Deformation 

Type:Total Deformati 

Unit:m 

Time:1 

2022/2/22 15:31 
1.9723e-7 max 
1.7532e-7 
1.534e-7 
1.3149e-7 
1.0957e-7 
8.766e-8 
6.574Se-8 
4.383e-8 
2.1915e-8 
0 min 


0.000 0.100 0.200 (m) 
ew 


0.050 0.150 


Dibox_shuzhi 

Total Deformation 

Type: Total Deformatié 

Unit:m 

Time:1 

2022/27122 15:35 
9.6994e-8 max 
8.6217e-8 
7,5439e-8 
60.4602e-8 
5.3885e-8 
4.3108e-8 
3.2331e-8 
2.1554e-8 
1.0777e-8 
0 min 


0.000 0.100 0.200 (m) 
B= = 


0.050 0.150 
7 齿轮 箱 位 移 云 图 


Fig.7 Deformation nephogram of gear box 
轴承 外 圈 故 障 频 率 对 轴 速 归 一 后 的 阶 次 由 式 

(20) 计 算得 到 。 
f BPFO 1 


f/f. 2 D, 
其 中 Z 是 轴承 滚动 体 数 目 。 


D, 
= z(1- ) =2.573% (20) 


通过 故障 轴承 处 加 速度 传 感 融 测 得 原始 信号 的 


频谱 如 图 8 所 示 。 


0.030 
0.025 
0.020 


0.015 


幅 值 /(m * s) 


0.010 


0.005 


0 i 10 15 20 25 
频率 /kHz 
图 8 外 圈 故 障 实测 信号 频谱 图 


Fig.8 Frequency domain spectrum of outer race 


fault measured signal 


图 8 中 在 0 ~1.8 kHz 的 范围 内 可 以 看 见 分 布 
较为 明显 的 共振 频带 ,低频 共振 频带 与 高 频 共振 频 


带 相 比 幅 值 较 大 。 为 验证 轴承 故障 信号 特点 ,分 别 选 


取 0 ~5 kHz,5 ~10 kHz 和 10 ~15 kHz 的 共振 频带 


行 带 通 滤波 求 包 络 谱 ,结果 如 图 9 所 示 。 


荫 进 


0.025 
区 
2 0020| Ty 0.021 029 8 
昌 0.015 
= 0010 X2.589 48 
旺 y 0.005 491 11 
是 0.005 
; 2 | 
阶 次 
(a) 0-5 kHz 共振 包 络 谱 
Wo 工 2.589 48 
wm 0.03 70.034 3699 
号 0.02 
旺 006 306 58 
至 0.01 .006 306 5 
0 2 四 6 8 rT 
阶 次 


fb) 5~10 kHz 共振 包 络 谱 


0 YN 
© 阶 次 
© (c) 10~15 kHz 共振 包 络 谱 
ca 图 9 外 圈 故 障 实 测 信号 包 络 谱 分 析 
© Fig.9 Envelope spectrum analysis of 
CD outer race fault measured signal 


Oz 过 对 比 发 现 ,0 ~ 5 kHz 的 包 络 谱 中 轴 速 对 应 
的 阶 次 占 主要 成 分 ,外 圈 故 障 对 应 的 阶 次 幅 值 较 小 ， 
无 明显 的 外 圈 故 障 特征 ,难以 通过 包 络 谱 判 断 轴承 
日 障 ;5 ~ 10 kHz 和 10 ~ 15 kHz 对 应 的 包 络 谱 均 
可 局 看 见 幅 值 较 大 的 外 圈 故 障 阶 次 ,外 圈 故 障 特征 
较 田 明显 。 与 10 ~ 15 kHz 共振 频带 对 应 的 包 络 谱 
相 区 ,5 ~ 10 kHz 共振 频带 对 应 的 包 络 谱 幅 值 更 大 ， 
特级 更 明显 。 这 是 因为 频谱 图 上 5 ~ 10 kHz 对 应 的 
共振 频带 幅 值 更 高 ,包含 的 冲击 信号 更 加 丰富 。 

对 比分 析 表 明 ,轴承 故障 冲击 会 激 起 系统 的 高 
频 共 振 。 对 高 频 共 振 频 带 的 包 络 信号 进行 频谱 分 析 
是 提取 轴承 故障 特征 的 有 效 手段 ,同时 也 说 明了 引 
入 谐振 器 模拟 系统 高 频 共振 的 重要 性 "1 。 

根据 实测 轴承 故障 信号 特点 ,设置 谐振 器 的 共 
振 频 率 为 7.5 kHz。 初 始 位 移 与 速度 均 为 0, 利 用 四 
阶 龙 格 库 塔 法 求解 动力 学 方程 组 得 到 各 个 自由 度 上 
的 位 移 、 速 度 和 加 速度 ,以 竖 直方 向 上 谐振 器 的 加 速 
度 作 为 响应 信号 进行 分 析 。 

仿真 计算 轴承 的 稳 态 响应 信号 加 入 信 噪 比 
15 db 的 高 斯 白 噪声 , 取 5 s 与 轴承 外 圈 故 障 实测 信 
号 归 一 化 后 进行 对 比 ,仿真 信号 与 实测 信号 的 时 域 
如 图 10 所 示 。 


ChinaX iv 合 作 期 刊 487 


幅 值 /(m * 8s) 


3 05 0 1 2 30 3 
时 间 /s 
(a) 仿真 信号 


幅 值 /(m * s”) 


2.5 
时 间 /s 
(b) 实测 信和 号 


图 10 ”外 圈 故 障 信号 时 域 图 
Fig. 10 Time domain spectrum of outer race fault signal 
从 时 域 信 号 图 中 可 以 看 出 ,仿真 与 实测 外 圈 故 
障 信号 在 时 域 上 都 具有 周期 性 冲击 。 实 测 的 齿轮 箱 
轴承 信号 大 致 上 与 仿真 信号 一 致 。 对 2.3 ~2.5s 的 
时 域 信号 进行 细节 对 比分 析 , 如 图 11 所 示 。 


X2.3734 
上 0.732 585 
* 


X2.353 83 
了 0.554 60 作 | 


pa HW 


3 A 15 ‘21 30 35 40 45 5.0 


上 2.393 18 
了 0.712 68 


幅 值 /m * s“) 


.0 
2.30 2.32 2.34 236 2.38 2.40 242 2.44 2.46 2.48 2.50 


时 间 /s 
(a) 仿真 信号 
08『 e235631 
0.597 902 
Y7397 11 
| y 0.436 63: 
XY2.3766 


了 0.408 808| 


幅 值 /(m *s) 


8 
2.30 2.32 2.34 236 238 2.40 2.42 2.44 2.46 2.48 2.50 
时 间 /s 
(b) 实测 信和 号 


图 11 外 圈 故 障 信号 细节 对 比 


Fig.11 Time domain signal detail comparison of outer race fault 
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6304 轴承 在 转 频 19. 66 Hz 下 的 外 圈 故 障 频 率 
在 时 域 上 对 应 的 理论 冲击 周期 为 0.019 8 s, 图 中 仿 
真 与 实测 信号 在 时 域 上 的 冲击 间隔 都 满足 这 一 特 
点 ,证 明了 仿真 信号 时 域 上 的 合理 性 。 

峭 度 (Kurtosis) 是 幅 值 归 一 化 的 四 阶 中 心 矩 ,也 
是 对 冲击 信号 敏感 的 统计 指标 。 谱 峭 度 ( Spectrum 
Kurtosis, SK) 能 够 反映 不 同 频率 下 瞬 态 冲击 的 强 弱 ， 
为 求 包 络 谱 中 心 频率 和 带宽 的 选择 提供 重要 参考 ， 
仿真 信号 的 快速 谱 峭 度 如 图 12 所 示 。 


在 4.6 级 别 上 /=13.933 5， 最 住 窗口 长 度 =48， 
中 心 固 率 =7.291 7 kHz， 带 宽 =2.083 3 kHz 


CO 
全 图 12 仿真 信号 快速 谱 峭 度 图 
Fig. 12 Fast-Kurtogram for the simulated signal 

“根据 仿真 信号 的 快速 谱 峭 度 图 上 的 最 佳 窗口 长 
度 绘制 信号 谱 峭 度 , 如 图 13 所 示 。 故 障 冲击 信号 在 
谐振 器 激 起 的 以 7.5 kHz 为 中 心 的 高 频 共振 频带 更 
加 是 富 ,证 明了 谐振 器 的 存在 能 够 模拟 故障 激 起 的 
系统 高 频 共振 ,与 实际 轴承 故障 机 理 更 加 吻合 。 

© 


峭 度 


0 5 10 15 20 25 
频率 /kHz 
图 13 仿真 信号 谱 峭 度 
Fig. 13 Spectrum kurtosis for the simulated signal 
为 进一步 验证 仿真 信号 的 准确 性 ,对 归 一 化 后 
的 仿真 与 实验 信号 进行 包 络 谱 分 析 , 仿 真 轴承 外 圈 
故障 信号 包 络 谱 和 实测 轴承 外 圈 故 障 包 络 谱 如 图 
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14 所 示 。 
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0.035 [1y 0.038 279 1| [0033 558 
XY7.723 16 


0.030 了 0.027 501 5 


全 0025 X10.2975 
70.021 6766 
三 0.020 X12.8719 
寺 70.016 8449 
至 0.015E IX1.577 19 
70.005 168 34 
0.010 X3.571 58 
yY 0.005 068 
0.005 
AAA 上 As 小 4 A 
0 2 4 6 8 10 2 4 6 
阶 次 
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阶 次 


(b) 实测 信和 号 


幅 值 /ms 


阶 次 


{c) 对 比 


图 14 外 圈 故 障 信号 包 络 谱 
Fig. 14 Envelope spectrum of outer race fault signal 
如 图 14(b) 所 示 , 实 测 信号 的 外 圈 故 障 特征 表 
现 为 包 络 谱 存 在 外 圈 故 障 特 征 阶 次 及 其 高 阶 谐 波 ， 
幅 值 随 阶 次 增 大 而 递减 , 且 由 于 轴 偏 心 距 的 影响 , 故 
障 特征 被 轴 速 调制 ,表现 为 存在 转 频 阶 次 间隔 的 劳 
辩 , 旁 兴 幅 值 相对 较 小 。 图 14(a) 中 仿真 信号 同样 


存在 对 应 的 故障 特征 阶 次 与 转 频 阶 次 的 旁 兴 。 

如 图 14(c) 所 示 , 与 实测 信号 相 比 ,仿真 信号 包 
络 谱 的 幅 值 较 大 ,主要 原因 是 传动 链 实验 平台 结构 
复 杂 ,存在 诸如 轴 偏 心 齿轮 路 合 等 耦合 振动 产生 的 
干扰 ;从 图 14(b) 的 包 络 谱 中 也 不 难 发 现 ,实测 信和 号 
存在 较 多 噪声 ;二 者 故障 特征 阶 次 值 吻 合 较 好 ,符合 
轴承 外 圈 故 障 的 振动 信号 特点 。 
由 式 (20) 计 算 可 得 外 圈 故 障 对 应 的 一 阶 特征 
阶 次 2.573 96 ,二 阶 特 征 阶 次 5. 141 92 三 阶 特征 阶 
次 7.721 88 四 阶 特征 阶 次 10. 2958 4 五 阶 特征 阶 次 
12. 869 8。 计 算 仿真 信号 和 实验 信号 的 特征 阶 次 与 
理论 值 的 相对 误差 如 表 2 所 示 。 

表 2 仿真 与 实测 信号 阶 次 相对 误差 


Tab.2 Relative error of order between simulated 


and measured signal 


相对 误差 /% 
仿真 信号 实测 信号 
0.0167 0.603 0 
0.016 5 0.603 2 
0.0166 0.6031 
0.0161 0.6028 
0.0163 0.6030 


水 , 而 实测 信号 特征 阶 次 与 理论 值 存在 一 定 的 偏 
差 呆 主 要 原因 是 轴承 在 实际 运行 中 ,滚动 体 在 滚 首 
中 并 非 做 纯 滚 动 运动 ,而 是 受 轴承 游 阶 和 润滑 条 件 
的 最 响 与 滚 道 存在 相对 背 动 。 

二 对 比分 析 发 现 ,仿真 信号 与 实测 信号 故障 特征 
相同 ,故障 特征 阶 次 相 一 致 , 幅 值 较为 接近 ,结果 表 
明 针对 齿轮 箱 轴 承 的 故障 建 模 是 准确 有 效 的 。 


3 结 论 


本 研究 基于 赫兹 接触 理论 ,分 析 深 动 体 与 故障 
的 局 部 接触 ,建立 了 集中 参数 的 轴承 故障 模型 ,引入 
谐振 器 以 模拟 轴承 故障 激 起 的 齿轮 箱 部 件 的 高 频 振 
动 。 对 实验 室 传动 链 平 台 的 齿轮 箱 进行 参数 识别 ， 
得 到 轴 和 人 齿轮 箱 在 水 平 及 竖 直 方向 上 的 刚度 。 以 
6304 轴承 外 圈 故 障 为 例 ,代入 各 项 参数 进行 仿真 计 
算 ,得 到 了 与 理论 轴承 故障 特征 相 一 致 的 结果 。 

将 实测 轴承 故障 信号 与 仿真 信号 进行 对 比分 
析 ,轴承 的 时 域 信号 及 包 络 谱 与 仿真 信号 吻合 较 好 。 
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计算 得 到 的 仿真 信号 故障 一 阶 特征 阶 次 与 理论 值 误 
差 约 为 0.016 7% ,证 明了 模型 具有 一 定 的 准确 性 和 
可 靠 性 。 该 模型 结构 较 简 单 且 计算 速度 快 ,针对 不 
同 尺 寸 的 轴承 能 够 快速 产生 相应 故障 信号 ,为 轴承 
故障 数据 库 的 建立 提供 了 便利 ,对 于 理解 轴承 故障 
冲击 振动 的 成 因 也 提供 了 分 析 参 考 。 
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